Radiation Losses of the Dominant Mode 
in Round Dielectric Waveguides 

By DIETRICH MARCUSE 

(Manuscript received March 5, 1970) 

The radiation loss theory that has been developed in a series of earlier 
papers is extended to the dominant mode of the round dielectric waveguide. 
The theory is applied to the calculation of radiation losses of abrupt steps, 
gradual tapers, and random wall perturbations of the round dielectric 
waveguide. 

The radiation losses caused by an abrupt step, and consequently the losses 
of tapers, are far higher for the dominant mode of the round dielectric wave- 
guide than they are for corresponding steps and tapers of the dielectric slab 
waveguide. However, the losses caused by infinitesimal random wall pertur- 
bations of the round waveguide are nearly equal to the random wall losses pre- 
dicted on the basis of the slab waveguide theory. In fact the losses of the 
dominate mode as well as the circidar electric TE 0X mode of the round rod due 
to random wall perturbations are very nearly the same. 

The theory is limited to circular symmetric distortions of the round dielec- 
tric rod {diameter changes). The radiation losses caused by steps of the 
round dielectric waveguide that carries the dominate guided mode have been 
verified by experiments at millimeter wave frequencies. 

I. INTRODUCTION 

A series of earlier papers was devoted to radiation losses of TE and 
TM modes in dielectric slab waveguides. 1-3 The radiation losses were 
assumed to be caused either by random perturbations of the waveguide 
boundary 1 or by steps and tapers of the slab waveguide. 3 Experiments 
to verify the radiation loss theory were conducted with millimeter 
waves in round teflon rods, and the theory was extended to cover this 
case. 2 

These earlier publications were limited to the simplified case of 
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electromagnetic fields that are independent of one coordinate. In the 
case of the slab waveguide we assumed 

T - ° (1) 

dy 

while 

i - ° (2) 

was required of the fields of the round dielectric waveguide. Restrictions 
(1) and (2) made it possible to separate the fields into transverse 
electric (TE) or transverse magnetic (TM) modes. 

The study of the simple slab waveguide yielded much useful informa- 
tion about the general properties of radiation losses and allowed us to 
infer the order of magnitude of the radiation losses caused by random 
wall imperfections. However, the dielectric slab is not a useful practical 
waveguide and can be used only as a simplified model to obtain informa- 
tion about the behavior of more realistic and more complicated struc- 
tures. Limitation (2) for the modes of the realistic and practical round 
dielectric waveguide precludes the application of the theory to the most 
important dominant mode of this structure. 

The present paper is devoted to a study of the radiation losses of the 
dominant mode of the round dielectric waveguide (optical fiber). To 
be able to handle the theory we still impose condition (2) on the deriva- 
tives related to the geometry of the waveguide but not on the field 
distribution. The resulting theory is still very complicated so that we 
must limit ourselves to sketching the theory and stating the final results. 

The radiation losses caused by random imperfections [obeying restric- 
tion (2)] are very nearly identical to the losses of the corresponding slab 
waveguide problem. However, the radiation losses of the dominant 
mode caused by steps and tapers in the waveguide are much higher 
than the corresponding losses of the TE or TM modes in the slab 
waveguide. The radiation losses of the dominant mode due to waveguide 
steps have been found experimentally to be in agreement with the 
theory. 

In order to allow the reader to obtain the information concerning the 
results of the theory unencumbered by complex mathematical formulas 
we start the paper with a discussion of the results. The remainder of 
the paper is devoted to an outline of the theory that was used to obtain 
these results. 
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II. NUMERICAL AND EXPERIMENTAL RESULTS 

2.1 Radiation Losses of Waveguide Steps 

We begin the discussion of the consequences of the radiation loss 
theory of the dominant mode of the round dielectric waveguide by 
considering the radiation losses caused by the abrupt step of the wave- 
guide diameter shown in Fig. 1. As described in Section II, the radiation 
losses caused by an abrupt step can be calculated by two different 
methods. The mode matching technique infers the loss from the trans- 
mission coefficient of the guided mode that continues to travel in the 
waveguide after it has passed the step. The radiation loss method 
accounts for the lost power by directly calculating the amount of power 
radiated into space. Both methods involve approximations so that we 
cannot expect to obtain exactly the same results either way. 

Figure 2 shows the results of both methods of calculation. The radia- 
tion loss caused by a step with a 2 /ai = 0.5 as a function of fca! (as 
computed by means of the mode matching technique) is shown as the 
dotted line in the figure, while the solid line represents the result of the 
radiation loss method. The curve holds for a dielectric rod with index 
of refraction n = 1.432 (n 2 = 2.05). This index was chosen since it is 
representative of teflon at a frequency of 55 GHz. The agreement of 
the two methods is remarkably good considering the approximations 
involved in deriving the theoretical expressions. 

Even better agreement is obtained by a similar calculation that applies 
to a dielectric rod with index of refraction n = 1.01 as shown in Fig. 3. 
Both figures are extended over kai values that correspond to single 
guided mode operation. There are other guided modes possible over 
part of the range of ka x values but these other modes do not couple 
to the dominant mode of the round dielectric rod because of the restric- 
tion on symmetry imposed by equation (2). It is in this sense that the 
operation of the waveguide is single mode. No other guided mode occurs 
under the imposed conditions. The shape of the two curves in Figs. 2 
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Fijf. 1- -Step in the round dielectric waveguide. 
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Fig. 2 — Relative radiation loss caused by an abrupt step with a 2 /ai = 0.5 of the 
waveguide. The two curves labeled dominant mode of the round waveguide were 
obtained by the mode matching technique (dotted line) and by the radiation loss 
technique (solid line). The two curves at the bottom of the figure labeled TE and 
TM modes represent the step losses of the slab waveguide. The radius a { (appearing 
in kai) belongs to the larger waveguide section. Index of refraction n = 1.432. 



and 3 is very similar. Both curves reach into high loss regions for small 
values of ka y . The curve of Fig. 3 is applicable to a clad optical fiber 
with 1 percent index difference between core and cladding. The curves 
shown on the bottom of Figs. 2 and 3 represent the step losses of TE 
and TM modes of the slab waveguide. 3 These curves are computed for 
the same index of refraction. The dimension a y (of kaO is the half width 
of the slab in the case of the slab waveguide. It is striking how much 
lower the radiation losses of the guided modes of the slab waveguide are 
compared to the dominant mode of the round dielectric rod. 

Because of the complexity of the theory and because the step loss 
results are so different for the round rod and the slab waveguide, it 
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appeared desirable to confirm the loss predictions of the theory with an 
experiment. The experiment was conducted with millimeter waves 
(approximately 55 GHz). A round teflon rod of 0.191 cm diameter was 
mounted between two metallic reflectors as shown in Fig. 4. The resulting 
resonant cavity could be excited through small holes in the reflector 
plates that, simultaneously, acted as supports for the teflon rod. Two 
teflon sleeves of 0.216 cm and 0.242 cm outer diameter could be slid over 
the teflon rods to produce a round dielectric waveguide with two steps. 
The losses caused by the steps could be determined from Q measure- 
ments of the cavity with and without the teflon sleeves. The results 
of these loss measurements (applied to one step) are shown as crosses 
in Fig. 5. This figure also shows the theoretical loss predictions of the 
mode matching (dotted line) and the radiation loss approach (solid line) 
of the theory. Note that the parameter value ka 2 = 1.1 of this figure 
uses the fixed value of the narrower portion of the waveguide as reference. 
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Fig. 3 — This curve is similar to Fig. 2 with n = 1.01 and a 2 /oi = 0.5. 



1670 THE BELL SYSTEM TECHNICAL JOURNAL, OCTOBER 1970 

| TO DETECTOR 



REFLECTOR 



,-TEFLON SLEEVE 



— TEFLON ROD 



REFLECTOR 



-TEFLON ROD 



J RECTANGULAR 
*~ WAVEGUIDE 



FROM KLYSTRON 



Fig. 4 — Experimental resonant cavity set up to measure radiation losses of wave- 
guide steps. 

The point a 2 /ax = 0.5 of Fig. 5 corresponds to the point ka x = 2.2 of 
Fig. 2. The measurements support the result of the round rod theory. 
The radiation losses of the slab waveguide even for much larger steps 
are still far lower than the measured values of these smaller steps of 
the round rod. 

It is not as easy to confirm the loss predictions of the slab theory 
since a dielectric slab waveguide is somewhat of an idealization. In 
particular it is hard to excite a slab with a mode that has no field varia- 
tion in the ?/-direction. In order to obtain some approximation to the 
slab waveguide we constructed a dielectric (teflon) ribbon whose 
dimensions on the narrower portion Avere 0.380 by 0.095 cm and whose 
wider dimensions were 0.380 by 0.190 cm. Note that only the narrow 
side is affected by the step. The losses of this ribbon waveguide with 
a 2:1 step were measured in the same resonant setup and compared to 
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the losses of a smooth ribbon with dimensions 0.380 by 0.095 cm. The 
radiation loss of the ribbon guide was AP/P = 0.08 for kd 2 = 1.1 
(or kdi = 2.2). This radiation loss value is shown as the circle in Fig. 5. 
It is apparent that the loss of the ribbon guide is far smaller than 
the loss of the round waveguide. It is about four times higher than 
the step loss predicted for the slab waveguide. However, we must keep 
in mind that the ribbon is only a poor approximation of the slab wave- 
guide. It is therefore not surprising that its radiation loss cannot be 
predicted by the slab waveguide theory. The slab waveguide apparently 
can tolerate steps in its width exceptionally well. 

2.2 Radiation Loss of Tapers 

The radiation loss theory that is presented in the theoretical part 
can be used to determine the loss of round dielectric waveguides with 
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Fig. 5 — Comparison of theory and experiment. The crosses are measured step 
losses of the round dielectric waveguide. The circle is the step loss of a ribbon guide. 
The carves represent the results of the mode matching theory (dotted line) and 
the radiation loss theory (solid line), (n = 1.432, ka* = 1.1.) Note that the curve 
parameter ha* uses the radius of the smaller waveguide Bection. 
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arbitrary diameter changes. Since the radiation losses of an abrupt step 
are very high for round dielectric waveguides it is interesting to study 
the radiation losses of gradual tapers. 8,7 

The calculation of the radiation losses of tapers can be simplified by 
observing that the dependence of /3 on the radius of the waveguide is 
nearly linear over a considerable range of values. Figure 6 shows the 
ratio of p /k as a function of ka for n - 1.432. It is apparent that a 
straight line approximation is possible in the region 1.2 < ka < 2.5. 

We study the radiation losses of two different tapers. The linear taper 
is the simplest and therefore the most reasonable taper to investigate. 
However, there are reasons to suspect that the linear taper may have 
higher radiation losses than other forms of tapers. It is apparent from 
equation (36) of Section II that the result of the integration (aside from 
the complicated factor I(p, z) which is difficult to evaluate) depends on 
the product of the derivative of the radius function a(z) with sine and 
cosine functions of the form cos f [0 o (z) - 0\dz. (j8„ is the propagation 
constant of the guided mode; is the z-component of the propagation 
constant of the radiation modes.) The oscillatory function has the 
tendency to cancel contributions from those functions that appear 
multiplied with it under the integrand. The more rapidly the cosine 
function oscillates, the more effective will be its canceling influence. 



1.28 

1.24 
1.20 




























































































1 no 

















0.2 0.6 



1.4 1.8 2.2 2.6 

ka 



Fig. 6 — Plot of the propagation constant /3 of the dominant mode of the round 
dielectric waveguide, (n = 1.432.) 
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This consideration shows that we would like to see the values of 
o (z) — as large as possible. The smallest possible value, and conse- 
quently the most harmful, is the value p (z) - k that is assumed at 
the upper end of the integration range in equation (34). However, 
because of the z dependence of p the values of /3 (z) — k are smaller 
at the narrow portion of the taper than they are on its wider portion. 
One might expect, therefore, that the narrow region of the linear taper 
contributes more to the overall radiation loss than its wider portions. 
It appears that the taper could be optimized if larger values of da/dz 
appeared at the wider end of the taper where the canceling effect of the 
sinusoidal functions is still more effective. Following this idea, it is 
possible to show that an exponential taper should distribute the radia- 
tion loss more evenly over its entire length in comparison with the 
linear taper. A linear taper and an exponential taper are shown in 
Fig. 7. The exponential taper was calculated from 

a(z) = a 2 + (fli - (h) exp ^ — 4.6 j 

This taper is designed to equalize the contribution of the integral (36), 
at least approximately, over the entire length of the taper assuming 
that I(p, z) is constant. The discontinuity of da/dz at z = does not 
contribute to the radiation loss. It would, therefore, be of no advantage 
to shape the taper such that da/dz is continuous over its entire length. 

The radiation losses of the linear and exponential tapers are compared 
in Fig. 8. Even though the radiation loss of the exponential taper is less 
than that of the linear taper, in agreement with our expectation, the 
amount of improvement is insufficient to warrant the greater com- 
plexity required to produce such a more complicated taper. Figure 8 
also shows that the radiation loss of a taper is far less than the losses 
caused by an abrupt step. The radiation losses can be made as small as 
desired with a taper of sufficient length. A linear taper with a length to 
waveguide radius (on the larger portion of the guide) ratio of L/a x = 400 
reduces the radiation losses, that would occur on an abrupt step, by a 
factor of 100. With X = 1 pm the value fca x = 2.5 is realized for o x = 
0.4 nm so that the taper would have an actual length of L = 160 pm 
or 0.16 mm. It is apparent that much longer, more effective tapers 
are feasible. 

Figure 8 indicates that there are two distinctly different regions. 
Below L/a l = 2 the taper is so short that it acts like an abrupt step. 
The beneficial effect of the taper makes itself felt only if the taper is 
long enough. The reduction of the radiation loss of a gradual taper 
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Fig. 7— The profile of the linear (dotted line) and the exponential (solid line) taper. 



compared to an abrupt step or steep taper is caused by the canceling 
influence of the (complex) exponential function in the integral of 
equation (36). 

2.3 Losses Caused by Random Wall Imperfections 

An important loss contribution is caused by the random deviations 
of the dielectric waveguide boundary from perfect straightness. These 
radiation losses have been investigated for the slab waveguide 1 and for 
the circular electric TE i mode. 2 The theory of radiation losses of the 
dominant mode of the round dielectric waveguide is sketched in 
Section III. 

We have seen that the radiation losses caused by arbitrary deforma- 
tions of the waveguide wall can be computed by describing the wall 
deviation as a series of infinitesimal steps. We have also seen that the 
single loss for large steps is far higher for the round dielectric waveguide 
than it is for the slab waveguide. We might thus worry that the losses 
caused by random wall perturbations may also be far higher for the 
dominant mode of the round dielectric waveguide. Fortunately, this 
pessimistic expectation is not true. The radiation losses caused by wall 
roughness of the round dielectric rod are no worse than they are for the 
modes of the slab waveguide. 

The random wall losses are treated on the basis of a statistical model. 
The correlation function describing the wall perturbation is assumed 
to be a simple exponential function that is characterized by two param- 
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eters, the rms deviation from perfect straightness A and the correlation 
length B. 

Figure 9 shows a series of curves of the normalized relative radiation 
loss as a function of the ratio of correlation length to waveguide radius 
B/a for a guide with index of refraction n = 1.432 (teflon). The curve 
parameter is the product of vacuum propagation constant times wave- 
guide radius, ka. Also shown for means of comparison is the loss of the 
circular electric mode of the round waveguide as a dotted line. It is 
apparent that the radiation losses of the dominant mode are approxi- 
mately equal to the radiation loss of the circular electric mode. A 
comparison with the results of Ref. 1 shows that the losses of Fig. 9 are 
approximately four times as high as the corresponding losses for the 
slab waveguide. For a meaningful comparison we must remember, 
however, that the slab waveguide losses were computed under the 
assumption that only one of the two slab boundaries was randomly 
perturbed. It seems reasonable to compare the losses of the round rod 
to a slab waveguide whose two walls are perturbed in a correlated way. 
In fact, if we assume that the thickness of the slab waveguide changes 
in a manner that provides equal but opposite displacement of each side 
of the guide we would obtain a four times higher loss than is shown in 
the curves of Ref. 1. The agreement between the radiation losses of the 
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j.'ig_ s — Relative radiation loss of the linear (dashed line) and the T exponential 
)lid line) taper, (n = 1.432, at/oi = 0.5, ka t = 2.5.) 
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Fig. 9 — Normalized radiation losses caused by random wall perturbations. The 
solid lines correspond to the dominant mode of the round guide, the dashed line 
represents the TEoi mode of this waveguide, (n = 1.432.) The curve parameters 
are the values of ka. 

slab waveguide and the random wall losses of the round dielectric 
waveguide is quite close. 

Figure 10 shows similar loss curves for a round waveguide with index 
of refraction n = 1.01. These curves too are about four times as high 
as the corresponding slab waveguide losses for the reason explained 
above. The curves of Fig. 10 are representative of the wall losses of 
a clad optical fiber with 1 percent index difference. As an example let 
us assume that we operate an optical fiber with a vacuum wavelength of 
X = 1 /jm. The value ka = 15 corresponds to a radius a = 2.4 /zm for 
the inner core of the fiber. If we assume that the correlation length of 
the exponential correlation function assumes its worst possible value 
B/a = 2.0, we find from Fig. 10 the normalized loss 



A loss factor of 



1 AP 
a = 7 =r = 2.3 km" 



= 10dB/km 



would be caused by an rms deviation of the waveguide radius = A 
9 10~ 8 cm = 9 A. This example shows how very stringent the 
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tolerance requirements can be. In a realistic case there will not only be 
variations of the waveguide radius. In addition we do not know the 
statistical model of the correlation function that must be applied in 
each case. However, comparison of different correlation function models 
has shown that the peak and its location in Figs. 9 and 10 is not de- 
pendent on the assumed statistical model. The decay of the loss curves 
toward increasing values of B/a is strongly model dependent. 



III. THEORY 

3.1 The Dominant Guided Mode 



The field components of an arbitrary guided mode in the waveguide 
are described by the following equations: 6 



E, = AJ v (i<r) cos vtf> 
H t = BJ v {nr) sin v§ 

E T = —\ K0 o AJi(i<r) + conB - J„(Kr) cos v<i> 
E* = \\ & A - J w (kt) + K(dnBJ' v (nr) sin v<f> 



(3a) 
(3b) 

(3c) 
(3d) 
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Fig. 10 — These curves are similar to Fig. 9 with n = 1.01. The curve parameters 
are the values of ka. 
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K 



Ha ~ 9 



n 2 ux A - J,(kt) + KfioBJifyr) 



rfmtoAJ'AKr) + \%B - J .(en) 



sin v<f> 



cos v<f> . 



(3e) 



(3f) 



These equations describe the field inside of the round dielectric rod, 
r ^ a. The functions J v are the Bessel functions of order v, a prime 
indicates the derivative with respect to the argument (not with respect 
to r). The parameter v must be an integer in order to make sine and 
cosine periodic functions of the aximuth <£ with period 2ir. The factor 

««"«-*" (4) 

was omitted from equations (3). The propagation constant /3 is related 
to the constants k and the free space propagation constant k by the 
relations 



k 2 = 



CO 6nMn 



and 



« a = n 2 k 2 - 0l , 



(5) 



(6) 



where n is the index of refraction of the dielectric material. The con- 
stants A and B are not independent of each other. Their mutual de- 
pendence is given by the boundary conditions for the field components. 
The fields on the outside of the dielectric rod r ^ a are given by the 
equations 

E z = CH?\iyr) cos v<f> 



H, = DH\ l \iyr) sin v<f> 



." Tilth 



E, - ■ -2 vyPoCH™' (iyr) + conD - HY\iyr) 
7 _ r . 



7 



o C-Hl l) (iyr) + iyapDH™'(vyr) 

r 



H r =~ 2 

7 



a>e C-Hl l) (iyr) + iyp DHl l) '(jyr) 

r 



H+ — — 

7 L 



iyueoCH^'iiyr) + p D -Hl l) (iyr) 



COS 1/0 

sin v<f> 
sin «£ 
cos «/> 



(7a) 
(7b) 

(7c) 
(7d) 

(7e) 

(70 



where ffj" is the Hankel function of order v and of the first kind. The 
prime indicates again its derivative with respect to its argument. The 
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argument is imaginary in order to ensure that the field distribution 
decays exponentially at large distance from the rod. The time and 
z-dependent factor (4) has again been suppressed. The parameter y is 
related to the propagation constant /3 by the equation 

y- = tf - lr- 

The field components were written down quite generally for an arbitrary 
guided mode. The lowest order or dominant mode of the guide follows 
from these equations with 

v = 1. (8) 

The following discussion will be limited to the special case v = 1. The 
connection between the amplitude coefficients and the determination of 
the propagation constant follows from the boundary conditions for the 
field components. The requirement that E, , E+ , H, and H+ are con- 
tinuous at the boundary r = a leads to the following eigenvalue equation 
for the determination of the propagation constant j8 of the guided mode 
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A few numerical values obtained from (9) are shown in Table I. The 



Table I — Some Numerical Values of n 
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connection between the amplitude coefficients as a consequence of the 
boundary conditions is stated in the following equations: 



(ka)(KaY 









(10) 



J x (ko) 



c = SFg^ < u > 



■7i(/ca) 
#"(*ya) 



g = r ,r, -B- (12) 



It is necessary to know the relation between the amphtude coefficients 
and the power P carried by the mode: 



5 fe {(or) W« + ■%»)] - 2J?(ica) 



-!- ■-'r-/;<M,)(i + ^5 



/io^B W+nV /3g + fc : 



+ 2 lti i r^ - ^;^jg) *■■ (l3) 

Equations (3) through (13) provide a complete description of the guided 
modes of symmetry cos 4>. The lowest order solution of the eigenvalue 
equation (9) is the dominant mode of the round dielectric rod. This 
mode does not experience a cutoff. In principle it can be supported by 
any round dielectric rod of arbitrarily small cross section and arbitrarily 
low frequency. All other modes of the round dielectric waveguide exist 
only above their respective cutoff frequencies. All entries in Table I 
belong to single mode (with cos <f> symmetry) operation. 

3.2 Radiation Modes of the Round Dielectric Rod 

The number of guided modes that the round dielectric rod can 
support is finite at any given frequency. In order to obtain a complete 
set of normal modes of the structure we need to consider also the 
continuous spectrum of unguided modes. 

Any solution of Maxwell's equations that satisfies the boundary 
condition is called a mode if its z-dependence (and time dependence) is 
given by equation (4). The guided modes are distinguished from the 
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unguided or radiation modes by the fact that their field distributions 
decay exponentially for increasing values of r outside of the waveguide. 
The radiation modes, on the other hand, extend to infinity. As their 
name indicates they are necessary to describe the radiation field outside 
(and inside) of the dielectric waveguide. Since there is no need to limit 
the functions describing the radiation modes to those that decay 
exponentially in the limit of large values of r we use a combination of 
Bessel and Neumann functions to express the unguided modes. How- 
ever, we must require that the field remains finite on axis at r = 0. 
These considerations allow us to express the unguided solutions of 
Maxwell's equations as follows: For r ^ a 

E, = FJ y (ar) cos «/> (14a) 

H t = GJ v {<rr) sin v<f> (14b) 

E r = -~ IrfiFJiW + "V-G ~ JM)\ cos v<}> (14c) 

E t = \ \pF V - J. (or) + ™nGJi(cr)\ sin v<t> (14d) 

H T = -\\ n 2 coe F - J, (or) + <rPGJ' p (or)\ sin v<j> (14e) 

H+ = -\ ln 2 aa>eoFJi(ar) + (3G V - J, (or) J cos v<p. (14f) 

There is now no restriction to the possible values that the propagation 
constant /3 can assume. The relation between and a is given by 

a 2 = n 2 k 2 - /3 2 . (15) 

The field outside of the dielectric rod, r ^ a, is given by 

E. = [HJ v (pr) + IN,(pr)] cos v4> (16a) 

H. = [KJ,(pr) + MN,{pr)] sin v<t> (16b) 

E T - -\ ipfi[HJ' r (pr) + INi(pr)] 

9 \ 

+ «m - [KJ,(pr) + MNXpt)]} cos »0 (16c) 

B t = -A(3- [HJ,(pr) + IN.(pr)] 
P I r 

+ po>n[KJl(pr) + MNi(pr)] \ sin v* (16d) 
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H r = -\ Le - [HJXpt) + IN,(pr)] 
p I r 



+ pftKJiipr) + MN'Apr)] > sin W> (16e) 



7 



P 



+ j8 V - [KJ,(pr) + MiV r (pr)]| cos v* (16f) 



with 



P 



= k 2 - /3 2 . (17) 

The Neumann functions N, are here expressed in the notation of 
Jahnke-Emde. 4 The determination of the coefficients of the radiation 
modes is complicated by an interesting phenomenon. The boundary 
conditions provide us with four equations. However, there are six 
undetermined coefficients in the set of equations (14) and (16). Even 
allowing for the fact that the power of the mode can be chosen arbitrarily 
so that one coefficient must remain undetermined by the boundary 
conditions, we have still one more coefficient than the boundary condi- 
tions, combined with the requirement of total power carried by the mode, 
are able to determine. This situation means physically that the sets 
of equations (14) and (16) represent a superposition of two modes that 
could be taken apart. A similar situation would have arisen in the case 
of the slab waveguide had we not been careful to separate the modes 
into even and odd field distributions from the very beginning. The 
present structure does not lend itself to a natural separation of the 
modes into even and odd ones. However, the formal field expressions 
(14) and (16) do, nevertheless, represent a superposition of two possible 
sets of modes. One might try to take arbitrarily either the coefficient 
F or G appearing in equation (14) equal to zero to try to separate out 
the two sets of modes. This procedure is mathematically beyond re- 
proach but it suffers from a practical inconvenience. The resulting sets 
of modes would not be orthogonal. It is very desirable to choose the 
modes in such a way that they are all mutually orthogonal to each other. 
It is therefore necessary to determine the coefficients in a way that 
assures the orthogonality of all the modes. The boundary conditions 
combined with the requirement of mode orthogonality and a certain 
amount of power carried by each mode are still not enough to assure 
a unique solution of our problem. This is not surprising since it is always 
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possible to combine two arbitrary vectors in an infinite number of ways 
into two mutually orthogonal vectors. 

The boundary conditions alone yield the following relations between 
the coefficients 



H = 1 (P«){[ 



2 P 



- n - J:(aa)NXpa) \F 
(n 2 - 1)A- 2 a v 



+ 



2 P 



P - J ,(aa)N r (pa)G 

pa coe a 






(18) 



I = % ( P a)< n 2 * Ji(aa)J,(pa) - J,(*a)J',(pa)\F 

(n 2 - l)k 



- J ,(ao)J Xpa)G 
pa coe u a 



(19) 



K = %{pa)\ ( £-^-[i{j .-r.Sipn-.r 



pa cop 



J r (<ra)N',(pa) - - J' p (aa)N p (pa) 



(20) 



M = | ( P a)\- 



(n* - l)fc s 



- J \(aa)J \,{pa)F 

pa cop CI 



+ [j J[{aa)J v (pa) - JA'aVlMjG]- (21) 

Equations (14), (16) and (18) through (21) are sufficient to satisfy 
Maxwell's equations and the boundary conditions. The coefficients F 
and G are, so far, completely arbitrary. We consider now two sets of 
radiation modes. The first set is distinguished by using the coefficients 
with subscripts F t and G'i while the coefficients of the second set are 
designated by F 2 and G- 2 . The two sets of coefficients must now be 
adjusted to render the two sets of modes orthogonal. One of the infinitely 
many solutions of this problem is 



El 
G. 






The ratio of F\/G\ is now no longer arbitrary but is given by 

F\ = faVr (g - b) 2 -h(e- df + (c 2 + f) 1» 
G x \J 1(9 ~ n 2 bf + (e - n 2 d) 2 + (c 2 + f)\ 



(22) 



(23a) 



p= = 
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with 

b = £ JUaaWtlpa) (23b) 

a 

c = & ~ V k §. jfaWfa) (23c) 

pa a 

d = R J' 1 ( f ra)J l (pa) (23d) 

a 

e = J^aajJiipa) (23e) 

= (n»- l)fc ^ j lM j l(pa) (23 f) 

per a 

e = J&aWKpa). (23g) 

Equation (23) was already specialized to the mode of symmetry cos <f>, 
taking v = I. The power carried by the radiation modes is given by 

■ Co 

The normalization of the radiation modes involves the delta function in 
the same way as it did in the case of the slab waveguides. 

3.3 Radiation Losses Caused by a Step 

It has been shown previously 3 that the radiation losses of arbitrary 
deformations of dielectric waveguides can be calculated from the 
knowledge of the radiation loss of a step. For simplicity we limit the 
discussion to waveguide imperfections that do not violate the condition 
(2). Condition (2) restricts the waveguide deformations to symmetrical 
changes of the waveguide diameter. More general deformations are far 
more difficult to calculate. 

A step in the round dielectric rod is shown in Fig. 1. We restrict 
ourselves to a dominant mode waveguide. The radius of the larger part 
of the waveguide must be small enough to ensure that only the dominant 
mode of the structure can propagate. Waveguides with larger radii 
suffer conversion losses to other guided modes in addition to the radiation 
losses. Such losses have been studied for the case of the slab waveguide 1 
and for circular electric modes in round dielectric waveguides. 2 

The radiation field can be expressed as an integral over all the radia- 



<&$*[> 



\c + (g-b)(*) S| -|/ + (e-d)^J f 



21 

F\ (24) 
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tion modes. Indicating the modes by script letters with the superscript i 
for the incident guided mode, r for the reflected guided and radiation 
modes and t for the transmitted guided and radiation modes we can 
write the boundary condition at the step as follows: 

sj° + a r &r + P fa r (p)6i: } (p) + p r (p)ej: , ( P )] d P 

Jo 

= c8!° + f [« ( (p)CG0 + p,(p)C(p)] dp (25) 

Jo 

s}« + a r sr + f [^gocm + p,(p)sj: ) (p)] dp 

Jo 

- d sj° + f Mp) &£( p ) + p,(p) Cto] dp (26) 

Jo 

3C<° + a r 3C< r) + f M [q r (p)Xl:\p) + pApW^p)] dp 

Jo 

- c t 3C;° + f [g ( (p)3C^(p) + p,(p)3Cj: ) (p)] dp (27) 

Jo 

3eJ° + <uC + f [gr(p)3ci:'(p) + Pr(p)3C( P )] dp 

Jo 

= c t 3ej n + f to,(p)3Ci: } (p) + p,(p)3Ci: } (p)] dp. (28) 

Jo 

These equations express the continuity of the transverse electric and 
magnetic field components at the step. The field components that are 
shown to be functions of p belong to radiation modes while field com- 
ponents that are not explicitly indicated as functions of p belong to the 
dominant guided mode. The amplitude of the incident guided mode is 
unity. The approximate solution of the equation system (25) through 
(28) follows the same reasoning that was presented for the case of the 
slab waveguide. 3 The coefficient c, can be calculated by using the ortho- 
gonality of the waveguide modes to the right of the step. The modes 
to the right of the step are not orthogonal to the modes to the left of the 
step because of the different waveguide size. It is thus not possible to 
separate the coefficients q r and p r (which, incidentally, belong to the 
two orthogonal sets of radiation modes) from the coefficient a r of the 
reflected guided mode. This problem makes it impossible to obtain an 
exact solution of the equation system. We neglect the reflected radiation 
modes when we calculate the coefficient c, . This approximation is 
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justified by the fact that for large steps the radiation favors the forward 
direction so that q r and p r can be assumed to be small. For very small 
steps where the ratio of forward to backward scattered power can be 
expected to be more nearly unity we need not worry about the coeffi- 
cients of the reflected radiation modes since the modes of the two guide 
sections become more nearly orthogonal to each other. 

The transmission and reflection coefficients can thus be determined 
approximately with the result 



and 



with 



c, = 



a r = 



2IJ 2 



h - J, 



h = o \~2 (fii-A. i — a>/ifii)((MiiA a — 2 B 2 ) ttu)/? _ \ 
2 (r/ 2 tl l {iy 2 a 2 ) 

-s + -2)Ji(Kia l )H^\i'Y 2 a i ) - -5 J \(/tiaJH{ l] \iy 2 a») 
5-1 (P1A1 — o} t iB 1 )(ni t }e A 2 — /3 2 5 2 )J 1 (/c 1 a 2 )J r ,(K 2 a 2 ) 

K X K 2 



(29) 



(30) 



(a>e n 2 p 1 A 1 A 2 + co M /5 2 -Si-B 2 ) 



I (2 2\ 

K x K 2 \Ki — K 2 ) 

■[KiJi(Kia 2 )J u (K 2 a 2 ) — k 2 J o{ki(i 2 )J \{t< 2 a 2 )] 
H (o}e IS 1 A 1 A 2 + upp 2 BiB 2 ) TwTTTT 2 : 

72 « 1 (*72«2) 



1 



_K, + 72 



ia 2 J ^a^H^ {iy 2 a 2 ) — *a,J r 1 (Kia 1 )//i(z7 2 a 1 ) 



+ — [a 2 / (Kia 2 )ff" ) («72a2) - fli-ZoO^^"'^^)]) 

2 X JipCiOi) 



+ 



7i-7?ffi 1, (n 1 ai) 



iffi ,) (t7ia 1 )fl P i I, (n»a,) 



.•72 CT <i>/- 



7i 



(31) 
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and with 



_\K 2 72/ 7a J 

7,7 2 "i W2«2^ 



+ 



f a A A L flf?^ ^ifeO J j N 



7,72(72 - 7?) V-"- 1 - ^ "*^"*' ff»'(t 7l a.) //i"(H2a,) 
[i7^| l, 07,«,)^i"M,) - *Y^" , (*Tifli)^i 1) (*yiai)] 



+ -(n 2 a>e 1) ^,l,,U. + a, M /3,# 1 E 2 ) 
Jj(ic 2 aa) 



+ 



1 



a 2 J n (K l a 2 )H\ ,) (iy 2 a 2 ) - OiJofaaJHS (tYiOi) 



+ ^ [ia 2 ./ 1 (K I a 2 )#J ,, (i72a 2 ) - i'a 1 J 1 (K 1 a 1 )//o 1, (t'7 2 ai)] 

72 



(32) 



The indices 1 and 2 attached to the coefficients and parameters indicate 
that the corresponding quantities belong to the waveguide to the left 
of the step (index 1) or to the right of the step (index 2). The coefficients 
A and B are the amplitude coefficients introduced in equations (3), (10) 
and (13). The factor P in equation (29) is the power carried by the 
incident guided mode. It was assumed that the power of all the modes is 
identical. The actual power carried by the mode is accounted for by the 
expansion coefficients a r , q r , p T , q t , p t , and c, . The power coefficients 
appearing in equations (13) and (29) are also identical. 

The theory of the dominant mode of the round dielectric waveguide 
is far more complex than the corresponding theory of the slab waveguide. 
This explains why the slab waveguide is so much more convenient to 
use for studying the general properties of radiation losses. 

The radiation loss caused by the step is obtained from 



f — 



c, | 2 - | a r | 2 . (33) 

However, the same radiation loss can also be obtained by accounting 
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for the power carried away in the radiation modes. We can therefore 
write also 

^= f (U| 2 +|p| 8 ) i f l dfi. (34) 

jr j-k p 

The subscripts r and t have been dropped from the expansion coefficients 
p and q. Both reflected and transmitted radiation modes are automati- 
cally included by extending the integration range from —A; to A; so that 
backward as well as forward traveling waves are included. The factor 
| /3 \/p appearing under the integration sign arose from converting the 
integration variable p to /3. 

The theory becomes much simpler when we limit the derivation of 
the p and q coefficients to small steps. It was shown in the work on slab 
waveguides 3 that arbitrary deformations of the waveguide wall can be 
treated as a succession of small steps. Even abrupt tapers can be de- 
scribed this way. In the limit of small step height Aa we can write 

Aa = ^ Az. (35) 

az 

The expansion coefficients q r and q t can approximately be obtained 
from equations (25) through (28) by a method that has been explained 
in some detail in Ref. 3. 

q(p) = [iM f z e-' !i *-"*"*. (36) 

The subscript r or t of q is no longer necessary since q r corresponds to 
negative values of /3 while q t corresponds to positive values of /3. The 
derivation of q has been simplified by expressing quantities pertaining 
to the waveguide to the right of the step in terms of the corresponding 
quantities for the waveguide to the left of the step. This approximation 
involves an expansion of the field quantities in Taylor series keeping 
only the first two terms of the expansion 

F(a 2 ) = F( ai ) + (*5J Aa. (37) 

The orthogonality of the modes belonging to the same section of wave- 
guide can be employed to eliminate many terms from the expressions. 
The resulting expressions for I(p, z) is far simpler than it would be had 
we considered a large step. We obtain 
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Up,*) = 



■ipyP 



./,M 



• (ft, + /3)7p(we„.-l — + o>pB —j 



y.f„(pa) + 'P^ttt^— JApa) 
y 2 + p J 

gj"(m) 

yi\oKpa) -r ' 
a 



- — J 1 (pa) 
IP 



yNo(pa) + ip ^7i>)V.1\ ATi(pa) 



ff, ll, («7«) 



+ (/c 2 + 0o0) 



7 T P 

if + B 0)j.W) 

3a da / 



-— tf,(pa) 

7P 



/ , ajif . n ai\, T . v I 



(38) 



The derivatives of the amplitude coefficients H , I, K, and M of equa- 
tions (18) through (21) are taken by keeping F and G constant. The 
reason for this prescription is the fact that the terms containing deriva- 
tives of F and G disappear from the equations because of mode ortho- 
gonality. 



dH irp [j a 2 - n 2 p 



- J u {aa) 



N (pa) - — 
pa 



\ #i(pa)l 

>a J 



[2 J 2pV 

+ — - pa + n\pa - — gl 

\_pa \ aa /. 



.2 P 



J.MA^pa) 
(n- - l)k 2 



+ [n 2 ^ - 1 J i (aa)N () (pa)>F + 



<rJ {oa)N x (pa) + pJ.MNoipa) - - J 1 (aa)N 



\(pa)^GJ (39) 



da 



= -^ [{a ** ~ wV / o (aB)[/.0ia) - j~ Q ^,(pa)] 
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+ 



B 



— pa 4- n I pa — — 
pa \ a<r 



J l (aa)J 1 (pa) 



- (n 2 * 5 - l)./ 1 M/„(pa)^ + ^ ^ 

coe po" 



aJ u (a-a)J,(pa) + p.J i {ad)J n {pa) — - ./,(<ra)./,(p' 



a)}(?] 



(40) 



dK 7rp , 2 -,w2 

da 2it 



| -^-VJoMiV^pa) + P .I l (aa)N l) (pa) - - ./,( ff a)Ar,(po)>/'' 
|_w^<rp I, o 



4- W (o-a) 



iV (pa) - - ^(pfl) 



+ — J.MN.ipa) - - ./, <ra)tf (pfl)>G 
paa a J . 

dil/ 7rp , 2 1 » 7 2 

— - = -— (n - l)/c 
da 2<r 



(41) 



jtifiap 

+ W (o-a) 



£- L/oWJ^pa) + p.J,{<ra)J {pa) - \ J 1 (<7a)J l (pa)>F 



■h(pa) - — Ji(pa) 
pa 



2 1 

+ JiiacfiJiipa) - - J x {<ja)J { (pa)\(j 

paa a 



(42) 



Equation (36) holds for q as well as for p. It is only necessary to insert 
Fx and O t in equations (38) through (42) to obtain the q coefficients 
while the p coefficients are obtained by replacing F t , 6?i with F 2 , G 2 . 
In order to use equation (34) for the relative power loss caused by 
radiation, it is necessarjr to calculate q and p with the help of equations 
(36) and (38). The coefficients appearing in these equations must be 
obtained from equations (39) through (42), and (10), (13), (22), (23), 
and (24). It should be apparent that this theory is of considerable 
complexity and can be handled only on an electronic computer. It is 
sad that the dominant mode in a round dielectric waveguide leads to 
such a complicated theory in comparison with the simple treatment of 
the slab waveguide. 
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3.4 Random Wall Perturbations 

An important source of loss is the radiation that is caused by small 
random perturbations of the waveguide wall. Such radiation losses have 
been discussed for slab waveguides in Ref. 1 and for round dielectric 
waveguides operating with the circular electric guided mode in Ref. 2. 
Equation (36) of our present analysis can be used to calculate the loss 
of the dominant mode of the round waveguide caused by random wall 
perturbations. Since the step losses of the dominant mode of the round 
waveguide are so much higher than the corresponding losses of TE and 
TM modes of the slab waveguide one might fear that the losses caused 
by infinitesimal random perturbations of the waveguide wall may also 
be substantially higher. Fortunately, this is not the case. The losses 
caused by random wall perturbations are of the same order of magnitude 
for all types of dielectric waveguides that have been studied so far. 

The losses caused by random wall perturbations are calculated with 
the help of a statistical model. Instead of using equation (34) for a 
particular waveguide we form the ensemble average (AP/P) over many 
statistically similar systems. For very slight perturbations of the wave- 
guide wall we can assume that I(p, z) is independent of the z coordinate 
and write equation (36), after a partial integration, in the form 

q(p) = +i(0 o - 0)Kp) f a(z)e- i ^~' )z dz. (43) 

•'o 

The argument z has been dropped from I(p) since this function is no 
longer dependent on z. The partial integration had the beneficial effect 
of causing a(z) instead of its derivative to appear under the integration 
sign. It was shown in Ref. 1 how substitution of equation (43) in (34) 
makes the scattering loss dependent on the correlation function 

R(u) = (a(z)a(z - a)) (44) 

after the expectation value has been taken. It is, therefore, possible to 
write the average value of the relative radiation loss as 

(ff) = 2L £ (fi - OT| I CI) G0 I 2 + I I C2, (p) IW) -^ dp (45) 
with 

F(0) = f R(u) cos 03 o - fiudu. (46) 

•'0 

The superscripts 1 and 2 indicate that the function I(p) has been com- 
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puted for both types of radiation modes that are associated with F t , G\ 
and F 2 , G 2 . 
If we use for the correlation function a simple exponential function 

R(u) = A 2 exp (-^) , (47) 



F{0) specializes to 1 

A' 



Fifi) = -r 
B 



(ft, ~ /3) 2 + ~ 2 



(48) 



IV. CONCLUSION 



We have found that the radiation losses of the dominant mode of a 
round dielectric waveguide are much higher than the corresponding 
losses of TE and TM modes of the slab waveguide. The radiation losses 
of the dominant mode of the round dielectric waveguide with an abrupt 
step have been verified by a millimeter wave experiment. The step 
losses of a ribbon waveguide were also measured and found to lie between 
the losses of the dominant mode of the round waveguide and the TE 
mode losses of the slab waveguide, but closer to the latter. It is thus 
apparent that the slab waveguide can tolerate abrupt steps exceptionally 
well. 

The radiation loss of a tapered round waveguide can be minimized 
by using a gentle taper instead of an abrupt step to accomplish the 
change of the waveguide radius. The losses of a linear taper are only 
slightly higher than the losses of a taper that was designed to equalize 
the loss contributions from different parts of the taper. It appears, 
therefore, that the design of optimum tapers is not profitable compared 
to their greater mechanical complexity. 

The losses caused by slight random wall imperfections are very similar 
for the dominant mode and the circular electric TE i mode of the round 
dielectric rod as well as the TE and TM modes of the dielectric slab 
waveguide. This result is surprising since the step losses of the dominant 
mode of the round waveguide are so much higher than the step losses 
of the slab waveguide. However, this result shows that the radiation 
losses caused by slight random wall perturbations can be studied with 
the help of the simple model of the slab waveguide and the results so 
obtained can be used to evaluate the performance of round dielectric 
waveguides. 
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